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2^ '■ Abstract 

^ I We investigate the connection between relativistic potential models for 

quark-antiquark bound states and the nonrelativistic models that have been 
used successfully to fit and predict the spectra of relativistic systems, as 
in the work of Martin. We use Martin's operator inequality \/p'^ + < 
(p^ + + m^)/2M to motivate the approximation of the relativistic kinetic 
energy terms in the spinless Salpeter equation by expressions of the nonrel- 
I ativistic form M + e + /2M for each quark. To investigate the validity of 

\^ ' the resulting approximation numerically, we generate energy spectra for qq 

■ mesons composed of two light or two heavy quarks using the spinless Salpeter 

. equation with the linear-plus-Coulomb potential typical of phenomenological 

fits to qq data, and then fit the lowest few states of each type using the ef- 
QhI fective Schrodinger description with the same potential. We find good fits to 

Q^' the lowest four calculated cc and the lowest three ss states either taking M 

^ ■ fixed at the value Mq = {p^) + that minimizes the Martin bound, or 

allowing Mq to vary in the fit. The energies of the lowest few cs states are 
then predicted with similar accuracy. The reasons for the success of the non- 
^ \ relativistic approximation are identified, and explain the success of Martin's 

' nonrelativistic predictions for the spectra of relativistic light-heavy mesons. 

However, we note that the agreement between the nonrelativistic and relativis- 
tic wave functions is not good, a point of potential concern for the calculation 
of transition matrix elements. 
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I. INTRODUCTION 



The development of potential models to describe the energy spectra of mesonic and 
baryonic systems has proved extremely successful. Phenomenological models that use a 
simple relativistic kinetic energy term and a scalar potential that incorporates the linear 
confinement and the short-distance color-Coulomb interaction suggested by QCD give good 
descriptions of the observed spectra of both heavy- and light-quark mesons and baryons |]T]-§| • 
Moreover, Duncan, Eichten, and Thacker have demonstrated a nontrivial connection 
between the relativistic potential models and rigorous numerical results from lattice QCD, 
showing that both the spectrum and the lattice wave functions for light-quark mesons are 
reproduced very well when the lattice potential is used in the relativistic wave equation 



+ mi + \l + m2 + V{r) 



^(r) = E^(r). (1) 



This equation, the spinless Salpeter equation, can be derived as a limit of the full Salpeter 
equation in which the "small-small" components of the Salpeter wave function are neglected 
and spin effects are averaged out as discussed, for example, in [Q. In the expression above 
p is the momentum of either quark in the center-of-momentum frame, mi and m2 are the 
quark masses, and V{r) is the effective potential between the quarks. 

We will be concerned here with the description of mesonic systems described as quark- 
antiquark bound states gQ, where the quarks q and Q may be the same or different. We as- 
sume that these systems can be described by the spinless Salpeter equation as demonstrated 
in [^, and will take V{r) as the linear-plus-Coulomb potential used in much phenomenolog- 
ical work. This also gives a good approximation to the lattice potential. For heavy quarks, 
the kinetic terms in Eq. (|I]) can be expanded in inverse powers of the quark mass to obtain 
the usual nonrelativistic Schrodinger Hamiltonian. This gives successful descriptions of the 
hb and cc states |^, even though the latter are close to being relativistic. More surprisingly, 
Martin showed that a nonrelativistic model based on a power-law potential could be ex- 
tended to include the clearly relativistic ss states, and was able using that model to predict 



successfully the masses of a number of then unmeasured light-heavy states ||10 . 

Although numerical methods have been developed which allow one to treat a relativistic 
kinetic term as easily as a nonrelativistic term ||3|JTT|-|r3| , it is important to understand why 
an ostensibly nonrelativistic treatment works and allows useful predictions to be made for 
relativistic systems as in the work of Martin and others. In this paper, we explore this 
problem theoretically, and develop a nonrelativistic approximation to the Hamiltonian in 
Eq. (0) based on an effective-mass expansion of the kinetic energy terms. We then study 
the accuracy of this approximation in reproducing the energy spectra and wave functions of 
relativistic qQ bound states by using the corresponding Schrodinger equation to fit "data" 
obtained by solving the spinless Salpeter equation. We find that it is possible to fit the energy 
spectra for the low-lying energy levels agree to within a few MeV for both heavy-heavy "cc" 
and light-light "ss" states. We are then able, using the nonrelativistic description, to predict 
the energies of the low-lying cs states to within 11 MeV. However, the effective quark mass 
M found in the fits is considerably larger than either the input quark mass m or the natural 
effective mass J {p"^) + m? expected from various arguments |ll^JTo| . We also obtain quite 



good fits to the relativistic cc and ss spectra, and good absolute predictions for the cs 



energies, using M = ^ {p'^) + rn? 
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We also study the wave functions in detail, and find qualitative agreement between the 
relativistic and nonrelativistic functions in the regions in which both are large provided the 
effective quark mass is used as a parameter in the fitting procedure. However, systematic 
differences are evident, and the nonrelativistic wave functions can be seriously in error locally, 
a problem that can limit the usefulness of the approximate wave functions in calculations of 
such quantities as transition matrix elements. 

In the next section, we develop the theory of the nonrelativistic approximation. We then 
outline the numerical techniques used to determine the energy spectra and wave functions, 
and discuss the results of the heavy and light fits and the light-heavy predictions in Sec. |T|, 
and summarize our conclusions in Sec. 0. 

II. THEORETICAL BACKGROUND 

Our objective is to approximate the relativistic potential model defined by the spinless 
Salpeter equation using a nonrelativistic Schrodinger description. Since only a small number 
of low-lying heavy- and light-quark bound states are actually known experimentally, an 
approximation will be successful for practical purposes if it reproduces the wave functions 
and the energy spectra for those limited sets of states. We will suppose that the potential 
V{r) is known and is kept fixed.Q This requires that our approximation for the kinetic 
energy terms in Eq. (|ID be accurate in some average sense for the low-lying states in both 
the heavy- and light-quark systems. 

Our nonrelativistic approximation to the kinetic terms is suggested by Martin's op- 
erator bound 

valid for an arbitrary mass M. The right hand side of this equation has the form of a 
nonrelativistic kinetic energy operator with an effective mass M, plus an additive constant 
that shifts the total energy. The equality in Eq. (|^) holds in momentum space at the 
momentum po = y/ M'^ — m^. Alternatively, the effective mass M is given in terms of the 
quark mass m and the point of tangency po of the curves defined by the two sides of the 
inequality by 

M^=m^+ pI (3) 

Because the Martin bound is an operator relation, the inequality in Eq. (|^) holds for 
expectation values in single states, and for averages of expectation values over sets of states. 
The choice p^ = {p^) would put the point of equality in Eq. (0) at the average value of p^ 
for the state or set of states under consideration. We would expect this choice for p^ to 



^Possible differences between the effective potentials V{r) for the heavy- and light-quark systems 
are outside our concern. However, we note that the potential is often varied in making phenomeno- 
logical fits to data on different systems, for example, in Q. This eases the problem of fitting the 
heavy- and light-quark systems together. 
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yield a reasonably accurate nonrelativistic approximation for the relativistic kinetic energy, 
a point noted in different contexts by other authors ||7;14;10]. More important theoretically, 



the effective mass M = y (p^) + m'^ minimizes the average value of the right hand side of 
Eq. (|^), so gives a least upper bound for the average of the relativistic kinetic energies when 
the average is calculated using the actual eigenfunctions for the relativistic problem. Using 
this value for M we obtain the relation 

p2 (^p2^ 



p^ + m^<M + - (4) 

The physical content of this result can be illustrated through a direct expansion of the 
square root operator. The standard expansion 

2 4 
P p 



p^ + m?' = m+ — - H (5) 

2m Sm-^ 

in powers of p^/m^ may be reliable for heavy-quark systems, but fails for light-quark systems. 
A possible solution to this problem is to consider an expansion about a fixed momentum p^, 



2M 8AP 



■ + ■■• (6) 



where M = \Jm? + p^. The expansion will give a good average approximation to the rel- 
ativistic kinetic energy provided the relevant values of p^ are concentrated near p^ with 
((p^ — Po)^) The numerator in this ratio has its minimum value for p^ = {p^). A 

comparison of Eqs. (0) and shows that the net effect of all the terms in Eq. (^) beyond 
the simple nonrelativistic result M + p^/2M is to decrease the kinetic energy. Note that the 
"relativistic correction" —{p^ — p^)"^ /8M^ to the kinetic energy operator in Eq. (^) does not 
have the standard form — p^/8M^, and would be expected to be much smaller in magnitude 
for pI close to (p^). 

To remove the strict inequality in Eq. in the following discussion, we will allow for 
an energy shift e' that includes the average contribution of the "relativistic corrections" in 
Eq. (^), taken as constant, and will use a nonrelativistic approximation to the relativistic 
kinetic energy operator of the form 



^Basdevant and Boukraa |14| consider the approximation obtained by including only the linear 
term in — (p^) in the expansion. A very different approximation which leads to a smaller effective 
mass M' = M/2 was proposed in and and studied in more detail by Lucha, Schoberl, and 
Moser This approximation was obtained by manipulating an inequality for matrix elements, 
Wp'^ + m"^) < \/Jp^y+TnF, and leads to an ambiguous result, in contrast to the operator inequality 
in Eq. (|2|). For example, the expression in Eq. holds as an operator inequality, but amounts to 
the addition of zero to the right hand side of the inequality for matrix elements when viewed at 
that level. The effective mass M' obtained in ||7|, p!^JT7[ ] is substantially too small, and the bound 



too weak, as will be seen in Sec. IIIB 
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p2 + m2^M + |^ + le, (7) 



where 

^ M M 4M3 



+ ■•■. (8) 



The content of this approximation is best illustrated in momentum space. In Figure |l|, 
we compare a model relativistic operator with m = 0.5 GeV, {p"^) = 3.75 GeV^, and M = 2 
GeV with the nonrelativistic approximation in Eq. (^. The curves corresponding to the 
relativistic and nonrelativistic expressions are tangent at = (p^). 

For all other momenta, the nonrelativistic approximation lies above the actual relativistic 
kinetic energy, as expected from the Martin bound. To improve the agreement between the 
operators for momenta away from the point of tangency, we can add a negative shift e' to 
the the nonrelativistic approximation as suggested above and shown in Figure Because of 
the negative curvature of the relativistic kinetic energy, it is also advantageous to increase 



the value of M relative to y (p^) + to move the point of tangency outward and reduce 
the slope of the nonrelativistic curve. This will be seen in our numerical results. The quality 
of the resulting approximation is evident in Fig. in which we compare the exact and 
approximate kinetic energies for the cc system over the region in which the wave function 
for the second excited cc state is large. The products of the kinetic energy operators with the 
squares of the momentum-space wave functions for the Salpeter and Schrodinger equations 



are compared in Fig. y. The details and interpretation of the fit are discussed in Sec. Ill B 



III. NUMERICAL INVESTIGATION OF THE NONRELATIVISTIC 

APPROXIMATION 

In this section, we will explore the accuracy of the nonrelativistic approximations derived 
above in the case of the cc, ss, and cs systems by comparing the results for the energy 
spectra and wave functions obtained by solving the corresponding Salpeter and Schrodinger 
equations. The relativistic kinetic energy operators will be approximated as in Eq. (|^, so 
that, for example, 

2 

= 2^p^ + ml + V{r) ^2M, + e, + ^ + V{r) (9) 
for charmonium. We will take a standard linear-plus- Coulomb form for V{r), 

V{r) = Ar-—, (10) 

with A = 0.203 GeV^ and B = 0.437. These values correspond to the potential parameters 
used by Fulcher for fits to the charmonium system [Q. We will concentrate on the L = 
states, and will consider the possibility of varying M as well as that of keeping M fixed 
at the value M = ^ {p^) + m? determined by a relativistic calculation. The best values of 
M and e in the Schrodinger equation, or of e alone, will be determined by making a least 
squares fit to the relativistic "data" calculated using the Salpeter equation. 
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A. Numerical methods 



We have calculated the relativistic energy spectra and wave functions using now-standard 
numerical methods developed elsewhere [jlT]-[T^. We first construct matrix representations 



for the potential V{r) and the positive operators Ef = + mf = — + in a suitable 
orthonormal basis of angular momentum eigenstates. The matrix Ef can be diagonalized 
by an orthogonal transformation U , Ef = f/Aj[/~^. The eigenvalues are necessarily positive. 

The square-root operator Ei = Jp^ + mf is then defined as UA]^'^U~^ where A^/^ is the 



diagonal matrix of the square roots of the eigenvalues [11,12|. With a finite basis, this 



construction reduces the solution of the Salpeter equation to the matrix eigenvalue problem 

{E^ + E^ + V- E) Ri = 0, (11) 

where Ri is the column-vector representation of the radial wave functions in the given basis 
for orbital angular momentum /. This equation can be solved by standard methods. 



As shown by Fulcher , the matrix elements needed in this construction can be calcu- 



lated analytically using basis wave functions 

^lm{^ = Riir)YUr) (12) 

with the angular dependence given by the spherical harmonics Yi^m and the radial wave 
function Rf{r) given by 

R^(r) = (3^/^{2(3rfe-^^'L^^+^ (2/?r) . (13) 
Here /3 is a length scale parameter and L'^'^'^ is the associated Laguerre polynomial [jl8 . 



This set has been investigated by several authors [p|,^pT|- p!3| , p!9|] . We find that a matrix size 
of 20 X 20 is sufficient to produce stable eigenvalues and wave functions. The same basis 
functions can be used to solve the Schrodinger equation as a matrix problem. 

In various figures which appear later, we will use the function Un,i{r) = rRn^i{r). The 
radial probability density for the quarks is just \un.i{r)\'^. We will also use the momentum- 
space wave functions (f)n,i{p) when analyzing our results. These are defined by the Fourier 
transform 

0n,i(p) = TT^ / drpjiprun,iir), (14) 
Zir^ Jo 

where ji is the standard spherical Bessel function. 

Finally, to determine the best values of M and e, we minimize the function 

N 

J2 {EB,k - E^R^kf (15) 

k=l 

for the lowest energy levels, varying M and e in the nonrelativistic Schrodinger equation 
with the calculated relativistic energies E^^j^ held fixed. 



6 



B. Results for heavy-quark systems 

We will use the cc system for our study of bound states of two heavy quarks. We use the 
quark mass rric = 1.320 GeV and the linear-plus-Coulomb potential determined by Fulcher 
jlSl in his Salpeter-equation fit to the charmonium spectrum. After calculating the exact 
Salpeter energy spectrum for those parameters to obtain our "data" , we fit the four lowest 
energy levels using a sequence of nonrelativistic approximations. Since it is frequently argued 
that the cc is almost nonrelativistic, we consider the standard Schrodinger kinetic energy 
2mc + p^/mc as well as the effective-mass approximation discussed above. In the latter case, 
we take either as fixed at the value {p^) + ml obtained using the Salpeter value of (p^) 
averaged over the states in question, or allow Mc to vary along with the energy shift ec- Our 
results are given in Table |. 

We see from Table | that the Schrodinger approximation is rather poor, with deviations 
of the fitted energies from the exact values ranging from 57 MeV in the ground state to 
173 MeV in the third excited state. The Schrodinger energies are all too high, and increase 
much too rapidly for the excited states, with a total change in the deviation of +116 MeV 
over the states considered. The failure of the Schrodinger approximation is not surprising 
given the rather large mean momentum in the Salpeter cc states, (jP') = 1.021 GeV^, where 

{p') = \j:{n\p'\n). (16) 

^ n=l 

This corresponds to a root-mean-square velocity (f^)^/^ = 0.61 for the quarks, and the 
system is semirelativistic. 

The energies obtained using the approximation in Eq. (0) with Mc = \J (p^) + ml are 
substantially better, with deviations ranging from -18 MeV for the ground state to +21 
MeV in the third excited state. Moreover, the approximate energies increase less rapidly 
than those for the Schrodinger approximation, with an excess increase of only 39 MeV 
relative to the Salpeter energies over the four states shown. The overall fit is good. The 
improvement in the mean energy is the result of including the energy shift ec- The fiattening 
of the deviations is the result of the larger value of the effective mass, with Mc = 1.662 GeV 
rather than the input mass mc = 1.320 GeV. The fitted value of the energy shift, e = —669 
MeV, is close to, and smaller in magnitude than the average kinetic term —(jP')/Mc = —614 
MeV as expected from Eq. (|^). The extra shift is associated with the terms omitted in 
Eq. (D. 

Finally, if we allow Mc to vary along with e in the fitting procedure, we obtain an 
excellent fit to the relativistic spectrum, with errors less than 3 MeV and a root-mean- 
squared (rms) deviation of 2.12 MeV as shown in Table |I[ However, Mc is now quite large, 
Mc = 1.861 GeV, while ec = —1.009 GeV. The large value of Mc is needed to slow the growth 
of the nonrelativistic kinetic energy with increasing p, and improve its agreement with the 
relativistic kinetic energy as remarked earlier. However, the resulting effective mass is not 
directly related to the charm-quark mass m^. 

As shown in Fig. |^, the variable-M nonrelativistic approximation leads to a seemingly 
excellent result for the kinetic-energy density. However, the relativistic and nonrelativistic 
wave functions do not agree precisely even for this fit as seen either in momentum space 
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in Fig. ^, or in position space in Fig. ^ Some quantities of interest such as leptonic ||20 



and electromagnetic transition rates are sensitive to these differences, and the nonrelativistic 
model must therefore be used with care. 

The increase in the heights of successive peaks in the nonrelativistic position-space wave 
function relative to the relativistic wave function, can be understood on the basis of the 
relativistic WKB approximation |^ .0 In particular, the velocity of a nonrelativistic particle 
is larger semiclassically than that of a relativistic particle in the region near the origin where 
the color-Coulomb potential is large, so the particle spends less time in that region and its 
wave function is consequently smaller. Correspondingly, its wave function is larger near the 
outer turning point. 

In Fig. P we show the effect of varying the mass on the wave function for the second 
excited state of the cc system. The lower masses shown bracket the input value of mc, while 
the highest mass is close to that obtained in the variable- mass fit, M^. = 1.86 GeV. It is clear 
from the figure that the wave functions are quite inaccurate for the lower masses, are not 
especially good even for the large effective mass \J (p^) + ml = 1.66 GeV, or the mass 1.86 
GeV obtained in the variable-mass fit. The trends in the wave functions discussed above 
are also clearly evident. 

Finally, in Figs. |^ and ^ we compare the total energy densities u*Hu and the potential 
energy densities u*Vu for the second excited states for the Salpeter equation and the optimal 
nonrelativistic approximation with Mc = 1.861 GeV. The difference between the potential 
energy densities results entirely from the difference in the wave functions. The systematic 
difference between the wave functions shows up clearly in Fig. |^. 



C. Results for light-quark systems 

Our results for ss system of two light quarks are given in Table y. We have used a 
strange-quark mass = 364 MeV in these calculations following Fulcher [Q, but have not 
changed the potential as he did, preferring to keep the same potential as for the heavy-quark 
system so as to be able to treat both systems simultaneously and predict the cs spectrum. 
The results are actually rather insensitive to because (p^)^/^ 7AA MeV ^ m^. The 
system is clearly relativistic, with an rms velocity (f^)^/^ = 0.90 for the quarks. 

The energies obtained with the effective mass Mg = ^ {p^) + ml are reasonably good on 
the average, but the approximate energies again increase too fast relative to the Salpeter 
spectrum. The fit obtained when Mg is allowed to vary is excellent, with the energies differing 
from the Salpeter energies by less than 4 MeV for the three lowest states considered. The 
fitted value of Mg has essentially no relation to the input mass m^. 

Unfortunately, the wave functions obtained in this case are poor even for the best fit to 
the spectrum. We compare the Salpeter and approximate energy densities in Fig. ^. The 
differences are due mainly to differences in the wave functions. Even the kinetic energy 
densities shows significant pointwise disagreement in this case. 



Numerical calculations show that the approximation is rather good in this case. 
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D. Predictions for the light-heavy system 



We consider finally the light-heavy system corresponding to the relativistic Hamiltonian 
of Eq. (|l]) with mi = rric and m2 = mg corresponding to the masses used in the discussion 
above. We use the nonrelativistic Hamiltonian 

= + M. + ^ + ^ + i + + V(r) (17) 

obtained by replacing the square- root operators in Eq. (|T]) by the approximation in Eq. (^. 
The kinetic term is of the standard Schrodinger form with a reduced mass M = M^Mc/ {Ms + 
Mc) given in terms of the effective masses rather than the quark masses. For the purpose 
of making predictions, we will keep the energy shifts and the effective masses fixed at 
the values determined separately for the heavy- and light-quark systems. These quantities 
would all be expected to change somewhat in the light-heavy system. For example, the 
masses Mj = (p^) + mf that minimize the Martin bound on the total kinetic energy 
change because of the different value of (p^) in the light-heavy system. The value of this 
quantity averaged over the three lowest states is {p'^)cs = 835 GeV, a value intermediate 
between the values (p^)cc = 1-021 GeV and {p'^)ss = 0.744 GeV obtained for the heavy- and 
light-quark systems. The energy shifts are given to leading approximation by ~ ~{p^) /M^, 
so also change. However, the conditions for minimizing the bound make the kinetic energy 
stationary with respect to the masses and Mg. As a result, by the Feynman-Hellman 
theorem there is no first-order change in the energies for small changes in (p^). More 
physically, the original nonrelativistic approximations for the kinetic energy operators are 
already good over a wide range of momenta as shown in Fig. ^, so the effect of the changes 
on the spectrum is not expected to be large. 

Our predictions for the Salpeter energy spectrum for the light-heavy system are shown 
in Table |II|. If we use the fixed values of the masses, the energies of the four lowest cs states 
are predicted to within 36 MeV as shown in the table. We note that the ground state is 
predicted to lie at too low an energy as a result of the large negative value of the energy shift 
defined above. However, an examination of Tables | and |l| shows that the predicted ground- 
state energies of the cc and ss systems are also too small. The usual fitting procedure adjusts 
the energy shift to minimize the deviations between the theory and the input data over the 
set of states considered. If we consider instead adjusting the energy shifts ec and eg to fit 
the cc and ss ground-state energies exactly, a reasonable procedure phenomenologically, we 
predict the normalized energies given in the third row in Table |IT1[ The ground state is now 
predicted correctly. However, the energies of the excited states cs increase too rapidly. This 
too rapid increase was also present for the cc and ss states. We note in this connection that 
the energies of the cs states are very close to the average of the energies of the corresponding 
cc and ss states. 

If we use instead of the fixed masses the fitted values of the masses and energy shifts for 
the heavy- and light-quark systems, we predict the energies of the lowest three cs states to 
within 11 MeV as shown in Table |T|. The largest difference occurs for the second excited 
state. The fits to the cc and ss energies are already excellent, and there is no reason in this 
case to renormalize the energy shifts. The closeness of the predictions to the actual energies 
would be expected given the results obtained for the cc and ss systems. In particular, the 
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nonrelativistic approximations to the kinetic energy operators are good in the regions in 
which the momentum-space wave functions are large. However, the final position-space cs 
wave functions are again not accurate. 

IV. CONCLUSIONS 

We find that the apparent success of nonrelativistic models for relativistic systems can be 
understood in terms of an approximation to the relativistic kinetic energy operator motivated 
by the Martin bound in Eq. (|]). Although the physical content of the approximation 
can be understood in terms of an expansion of the relativistic operator about a mean mo- 
mentum squared Pq, given optimally from the bound as = (p^), the series expansion is 
not necessary. What is important is to obtain a good average representation of the kinetic 
energy operator of Schrodinger form. We observe in this connection that the approximation 
can be improved significantly by allowing an extra energy shift to eliminate the inequality, 
and, if desired, also allowing the effective mass M appears to vary. 

We have investigated the effectiveness of this procedure in detail by using the nonrel- 
ativistic approximation to fit "data" obtained by solving the relativistic Salpeter equation 
for the linear-plus- Coulomb potential used by Fulcher [|T^] in fits to the the charmonium 
spectrum. We find that the nonrelativistic approximation for the kinetic energy operator 
in Eq. (|^) gives generally good descriptions of the Salpeter energy spectra for the cc and 
ss systems, taken as examples of bound states of heavy and light quark pairs. The results 
obtained with the effective masses fixed at the values (p^) + suggested by minimizing 
the Martin bound over a set of states are good, but the excited state energies generally 
increase too rapidly if the potential is kept fixed. The results obtained when M is allowed 
to vary in the fitting procedure are accurate to a few MeV is all striking result. 

We believe that the theoretical understanding of the success of the nonrelativistic 
effective- mass approximation developed here provides a justification for Martin's nonrel- 
ativistic treatment of heavy- and light-quark systems, and explains the unexpected success 
of his predictions for the masses of light-heavy systems [P,^p!0| . 
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FIG. 1. We show the relation between the relativistic kinetic energy operator and the nonrel- 
ativistic approximation given by the Martin bound, Eq. (ph. The effective mass M and the quark 



mass m are related by 



m 



+ (jP). The values used are M = 2 GeV and m = 0.5 GeV 



which give the local equality at p"^ = 3.75 GeV^. We also plot the nonrelativistic approximation in 
Eq. (0) with an energy shift e/2 = -1.1 GeV instead of the shift -{p^)/2M = -0.94 GeV in the 
Martin bound. The agreement between the two expressions is improved at low and high momenta. 




0.25 0.5 0.75 1 1.25 1.5 1.75 2 2.25 2.5 
p (GeV) 

FIG. 2. We plot the relativistic and the approximate nonrelativistic kinetic energy operators 
with the square of the second excited state wave function for the cc system of Sec. IIIB| superposed 
to show the approximate agreement of those operators in the region in which the wave function is 
large. M = 1861 MeV, e = -1009 MeV, and = 1320 MeV. 
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FIG. 3. We compare the relativistic momentum-space kinetic energy density for the second 



excited state of the cc system of Sec. IIIB with the density obtained using the nonrelativistic 
approximation in Eq. (0) with M = 1861 MeV and e = -1009 MeV. 
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FIG. 4. We compare the Salpeter momentum-space probabiUty densities for the ground 

state and second excited state of the cc system of Sec. [II B| with the densities obtained using the 
nonrelativistic approximation for the kinetic energy given in Eq. (0) with M = 1861 MeV and 
e = -1009 MeV. 
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FIG. 5. We compare the Salpeter position-space probability densities for the ground 

state and second excited state of the cc system of Sec. IIIB| with the densities obtained using 
nonrelativistic approximation for the kinetic energy given in Eq. (0) with M = 1861 MeV and 
e = -1009 MeV. 
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FIG. 6. We show the effect of varying the mass M in the nonrelativistic approximation for 
the kinetic energy operator, Eq. (^), on the quark radial probability density |n(r)p for the second 
excited state of the cc system of Sec. |IIIB . The best agreement of the wave functions is achieved 
for the large effective mass M « 1.8 needed in fitting the energy spectrum. 
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FIG. 7. Plot of the radial energy density u*Hu for the second excited state of the cc system. 
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FIG. 8. Plot of the potential energy density u*Vu for the second excited state of the cc system. 
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FIG. 9. Plot of the energy density u*Hu for the second excited state of the ss system. The 
differences between the relativistic and nonrelativistic cases result mostly from differences between 
the wave functions. 
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TABLE I. Comparison of the exact Salpeter energy spectrum for the "cc" system of two heavy 
quarks with the spectra obtained in various nonrelativistic approximations. The Schrodinger ap- 
proximation involves the kinetic energy p^/mc- The nonrelativistic (NR) approximation is defined 
in Eq. (|^, and is considered both with the effective mass Mc fixed at the value ^^(p^y^Fm^, and 
with Mc allowed to vary, ec is the energy shift defined in Eq. (^). 

Model Mc(MeV) ec(MeV) Ei (MeV) E2 (MeV) E^, (MeV) E^ (MeV) l^E rms (MeV) 



Salpeter 


1320 




3067 


3668 


4112 


4486 




Schrodinger 


1320 





3114 


3755 


4241 


4659 


119 


NR, M fixed 


1662 


-669 


3049 


3660 


4116 


4507 


14.5 


NR , M free 


1861 


-1009 


3069 


3667 


4109 


4488 


2.12 



TABLE IL Comparison of the exact Salpeter energy spectrum for the "ss" system of two light 
quarks with the spectra obtained using the nonrelativistic (NR) approximation defined in Eq. (|^, 
taken either with the effective mass Ms fixed at the value V {p^) + rni or allowed to vary, is the 
energy shift defined in Eq. (P). 

Model Afs(MeV) e^(MeV) Ei (MeV) E2 (MeV) £^3 (MeV) /^E rms (MeV) 

Salpeter 364 — 1531 2222 2744 — 

NR, M fixed 828 -795 1503 2219 2775 24.2 

NR, M free 989 -1022 1533 2218 2746 2.83 



TABLE in. Comparison of the exact Salpeter energy spectrum of the heavy-light "cs" system 
with the spectrum obtained using the nonrelativistic approximation for the kinetic energy given 
in Eq. (p!7|). The normalized nonrelativistic spectrum with the masses Mj fixed is obtained by 
adjusting the energy shifts to match the ground states of the cc and ss systems exactly. 



Model 


El 


E2 




A£; rms (MeV 


Salpeter 


2319 


2957 


3438 




NR, M, e fixed 


2296 


2963 


3474 


24.9 


Normalized 


2319 


2986 


3497 


38.0 


NR, M free 


2319 


2961 


3449 


6.8 
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